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" God forbid that we should give out a dream of our own imagination 

for a Pattern of the World." 

— Francis Bacon, Novum Organum 



1 Introduction 

In this article we will be concerned with some possible physical realizations 
of quantum gates that are useful for quantum information processing. After 
a brief introduction into the subject, in Sec. 2 we will focus on a particu- 
lar way of using atoms in one-dimensional optical lattices as carriers of the 
quantum information. In Sec. 3 on the contrary, the information carriers will 
be photons that interact via effective nonlinearities which arises from mixing 
at passive linear optical elements and postselection through photodetection. 
These two seemingly different implementations have in common that their 
decoherence mechanism is described by a single theory, namely that of quan- 
tum electrodynamics in causal media which will be the subject of Sec. 4. 
Figure 1 should serve as an overview of the subject areas covered. 
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Fig. 1. Connections between the subject areas covered. 



2 Stefan Scheel, Jiannis Pachos, E.A. Hinds, and Peter L. Knight 

1.1 Why quantum information processing? 

Apart from the academically driven curiosity to learn more about informa- 
tion in a quantum-mechanical setting, there is very urgent practical need 
to investigate information processing and computing from a quantum the- 
ory point of view. Over the last 30 years the number of transistors on an 
integrated circuit, i.e. the complexity of the computer made up from those 
chips, doubles roughly every 18 months. This empirical behaviour is famously 
known as Moore's first law (after the co-founder of Intel Corp., G. Moore). 
An extrapolation reveals that by the year 2017 a bit of information will need 
to be encoded in a single atom. Even if Moore's law breaks down before we 
eventually come to this point, at the current growth rate, by the year 2012 
the dimensions of logical elements will be so small that quantum effects upon 
computation cannot be neglected any longer. 3 Thus, technological progress 
necessitates the study of what implications quantum theory has on compu- 
tation. 

Another reason for looking deeper into quantum computing lies in the po- 
tential ability to simulate the temporal evolution of (possibly chaotic) quan- 
tum systems. The intrinsic and potentially massive parallelism of quantum 
computers that stems from linearity of quantum mechanics and the result- 
ing superposition principle, would allow to investigate quantum Hamiltonian 
systems without the need for an exponential temporal overhead on a classical 
computer (see e.g. [1]). 

1.2 Quantum gates vs. classical gates 

"Information is physical." 
— Rolf Landauer 

The most important difference between a logical gate known from classical 
information processing and a quantum gate is its reversibility or conservation 
of information. In order to see what that means, let us consider for example 
the logical NAND gate whose action of logical 0's and l's are given by the 
truth table 1 . Equivalently, we can describe it by the action of X\ A X 2 on two 
Boolean variables X\ and X 2 - This gate is generic for classical gates in that 
it has two (or more) inputs and only a single output. In quantum mechanics, 
we are allowed to form linear combinations or superpositions of the logical 
basis states, i.e. co|00) + ci|01) + C2 1 10) + C3III). The formal application of 
the classical NAND on this quantum superposition would result in a state 
oc C3IO) + (co + ci + C2)| 1) ■ The resulting state contains far less information 
than the original state. Whereas before the gate operation all the weight 
coefficients Ci gave us information about the quantum state, after the action 

3 The most recent Intel Pentium 4 processor already contains features of the size 
of 90nm, i.e. 900A. The lithography used to generate such structures therefore uses 
XUV light! 
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Table 1. Truth table of the logical NAND gate 



inputs | output 


00 


1 


01 


1 


10 


1 


1 1 






of the classical NAND only the sum c + C\ + c 2 plays a role in determining 
the weight of the logical All the information about the mutual differences 
are lost during the operation. However, this state is obviously not properly 
normalized anymore. The formal action of a classical gate is therefore actually 
an ill-posed operation in quantum mechanics. In fact, there is no way of 
making sense of classical logical operations acting on quantum superpositions. 

In quantum information processing, all operations have to be unitary. 
This means on one hand that the number of output degrees of freedom must 
equal the number of input degrees of freedom, and on the other hand it 
means that no information is lost about an initial superposition. In general, 
we could allow even more general operations such as completely positive maps 
of which the unitary operations are a subset. However, not even using the 
correct number of inputs and outputs together with no loss of information 
is sufficient to define a quantum operation. For example, let us consider an 
operation in which a given quantum state co|0) + ci|l) is transformed into its 
orthogonal complement, hence into the state — cf |0) +Cq|1). This operation 
could be thought of as being the analogue of the classical NOT operation 
which is defined asInI:IM = 0. However, it turns out that the 
quantum-NOT operation is not a completely positive map and therefore can- 
not be implemented by any quantum circuit. In fact, the map describing this 
transformation is anti- unitary [2] with determinant —1. 4 

Universal set of quantum gates 

Now we have shown that there is no strict connection between classical logical 
gates and quantum mechanics, we will next briefly discuss which operations 
are compatible with quantum mechanics. Let us first look into operations 
acting on a single logical qubit state Co|0) + ci j 1) . The dynamical group 
associated with this state is the unitary group SU(2) whose action can be 
given in terms of its generators, the Pauli spin operators <7j (i — x,y,z), 
as U — e la (T . Here, a. denotes the direction of the rotation axis and the 
associated rotation angle. By decomposing U into exponential factors, one 
can derive the so-called Euler decomposition U = e Wl<T *e W2<T «e W3CTz which 

4 In this context, the notion of universal gate has been invented. A universal 
NOT-gate would be the one that, when averaging over all input state, comes closest 
to the NOT-operation. 
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just describes a sequence of elementary rotations. The Pauli operators have 
the following representation in the single-qubit basis: 

^ = |0}(1| + |1)(0|, (1) 
a w = i(|0)<l|-|l)<0|), (2) 
^ = |0)(0|-|1)(1|. (3) 

An equivalent set of operators would be a± — a x ±ia y and a z which results in 
an operator decomposition of the form U = e ibi(i+»z) e ib 2 » + ie b 3 a- ie b 4 (i-a z ) _ 
This result will be used later in Sec. 3 to describe beam splitters. 

By definition, the Pauli operators are enough to represent all other single- 
qubit operations. An important example is the Hadamard gate H (not to be 
confused with the Hamiltonian), 



(|0} + |1))(0| + (|0)-|1»(1| , (4) 



which transform logical states into coherent superpositions. From the defi- 
nition of the Pauli operators it is obvious that the Hadamard gate can be 
written as H = (a x + <t z )/\/2. 

As for two-qubit gates, the structure of the unitary operator associated 
with the dynamical symmetry is not so obvious. We will therefore restrict 
our attention to some particularly useful gates. One of the most important 
nontrivial two-qubit gates is the controllcd-phase gate C v whose truth ta- 
ble is given in Table 2. The controlled-phase gate (or its special case the 



Table 2. Truth table of the controlled-phase gate 



input output 



|00) 


|00) 


01) 


|01) 


10) 


|10) 


111) 


e*|ll) 



controlled-cTj gate with ip — tt), together with the set of Pauli operators, 
forms a universal set of quantum gates [3] by which we mean that all possible 
quantum networks can be build up from them. Therefore, if one is able to 
build these four gates, one can generate arbitrary quantum operations by 
concatenating them (Fig. 2). 

For example, another nontrivial two-qubit gate is the controllcd-NOT (or 
CNOT) under whose action the logical states 1 10) and |11) are exchanged and 
the other left untouched. 5 This operator can be realized by the combination 



5 The associated operator can be written as |0){0| <g> / + |1)(1| ® o z . 
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Fig. 2. A quantum network can be build from a universal set of quantum gates. 
Here we sketch a three bit register which is manipulated by a sequence of one and 
two qubit gates to realize a chosen unitary transformation which effects the relevant 
algorithm. 

HiC^Hi, i.e. by two Hadamard gates acting on mode 1 and a controlled- 
& z gate. Other two-qubit gates of importance are the swap in which the 
qubits are simply exchanged and the square root of swap which serves as an 
entangling gate. 

2 Atomic realisation — atom chips and the Mott 
transition in optical lattices 

There are many possible realizations of quantum bits: laser-cooled trapped 
ions, cold trapped neutral atoms, atoms in high-Q single-mode cavities as 
well as condensed matter candidates such as quantum dots or Josephson 
junctions. For any one of these to be viable, the effects of the environment 
(dissipation and temperature) must be minimized, if not entirely eradicated. 
This simply means we need to use cold state initialization. 

In this section, we will describe a possible way of implementing quantum 
computation with cold atom technology. This includes the application of 
optical lattices in a sufficiently cold cloud of atoms showing Bose-Einstein 
condensation (BEC). The lasers which generate the periodic spatially varying 
trapping potential through the AC Stark effect create an optical lattice strong 
enough to induce a quantum phase transition from the superfluid phase that 
characterises the BEC to the Mott insulator phase, characterised by a regular 
structure of one atom per lattice site that can serve as a quantum register. 
A universal set of quantum gates can then be realised by manipulations of 
the lattice potential with additional laser fields. 

2.1 Bose Einstein condensates and the Mott transition 

In recent years it has been realized that Bose-Einstein condensates (BECs 
for short) can undergo a phase transition if loaded into a three-dimensional 
periodic potential which for example can be realized by standing-wave op- 
tical fields [4]. That is, one starts off with a BEC in its superfluid phase in 
which the relative phases (or rather correlations) between the atoms are well- 
defined such that the whole ensemble of atoms can be described by a single 
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Fig. 3. Counter-propagating laser beams induce a periodic spatially varying trap- 
ping potential through the AC Stark shift. 

macroscopic wave function (in first approximation). By loading this conden- 
sate into the optical lattice (see Fig. 3) the number of atoms per lattice site 
is undetermined and can vary widely. The ground-state wave function for N 
atoms in a lattice with M sites is therefore 



where the a] denote creation operators of an atom at the lattice site i. How- 
ever, when increasing the strength of the potential by increasing the power 
of the laser beams that create the standing-wave potential, eventually there 
will be a phase-transition to a state of the condensate in which each lat- 
tice site is occupied by a fixed and well-defined number of atoms (ideally we 
would like to have exactly one atom per site) . In this so-called Mott-insulator 
phase the relative phases (or correlations) between neighboring lattice sites 
are undetermined. The ground-state wave function here looks essentially like 



where n is the number of atoms per lattice site. Experimental evidence of this 
phase-transition has been obtained in the beautiful experiments described in 
[5, 6]. Although a Bose-Einstein condensate really exists only in three dimen- 
sions (since only there we find a phase transition from a thermal cloud to a 
condensate) , there are analogous systems such as the quasi-condensate [7] and 
the Tonks-Girardeau gas [8] in one dimension that have similar properties. 

The effective interaction Hamiltonian that can be derived from the Gross- 
Pitaevskii equation under the assumption that the atoms have localized 
single-particle wave functions can be written as 




(5) 



M 



\*m) oc n(«i) n i°> 



(6) 



1=1 




-^n,(n 2 - 1) - J 



CH+i + oJ i+1 ai) . 



(7) 



This Hamiltonian is also known in the literature as the Bose-Hubbard Hamil- 
tonian. The first term with coupling strength U is the collisional energy of 
atoms occupying the same lattice site. Obviously, if there is zero or just 
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one atom per site, this term vanishes identically. The second term is the so- 
called hopping term which essentially is given by the overlap of the localized 
wavefunctions at neighboring sites and describes tunneling between adjacent 
lattice sites with tunneling strength J. If, for example, the optical lattice is 
formed by the standing wave of a one-dimensional cavity mode having width 
L, the trapping potential is given by 



V(x) — — Vq sin 2 kxexp 



2r2 
1? 



(8) 



where r denotes the transverse distance from the lattice axis and k the wave 
number of the standing wave. With this potential, the collisional strength 
can be approximated by 



U : 



4a„Vn E 



3/4 „l/4 



R 



XL 



(9) 



where a s is the scattering length of the atomic collisions and Er the atomic 
recoil energy. An analogous derivation shows that the tunneling rate can be 
expressed as 



J 



Er 



■ exp 





(10) 



It is now apparent that by changing the potential depth Vq one can tune 
the system of atoms into cither of the two (superfluid or Mott-insulator) 
phases as seen in Fig. 4. For example, if Vq is relatively small compared 



U/E R , J/E E 



V /E E 



Fig. 4. The values of U/ Er (solid line) and J / Er (dashed line) are plotted 
against Vq/Er for typical values of the relevant length scales appearing in Eq. (9) 
(o s = 5.6nm for 87 Rb, L = A = 10/im). 



to the recoil energy Er, the tunneling rate J will be large and the atoms 
will be delocalized and form a superfluid. Increasing Vq means exponentially 
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decreasing the tunneling rate and the atoms will become stuck in their re- 
spective potential wells and form the Mott-insulator phase. In Fig. 5 we show 




1234567B9 10 1 23456789 10 

Fig. 5. Distribution of 8 atoms among 10 lattice sites in the limits U/J <C 1 (left 
figure) and U/J S> 1 (right figure). 



an example calculation for the population distribution in a one-dimensional 
lattice with 10 sites occupied by 8 atoms. The figure on the left depicts a 
situation in which the collisional energy U is small compared to the tunnel- 
ing rate J (U/J <C 1) and hence a superfluid phase exists, witnessed by a 
Gaussian-like distribution. On the other hand, the figure on the right shows 
the one-by-one distribution of atoms for U/J^>~1 which closely resembles 
a Fock state at each site. In fact, the critical value for obtaining the phase- 
transition is U/J~ 11.6 [9]. Note that the calculation of the ground-state 
wave function can only be done numerically which amounts to computing 
the eigenvector corresponding to the smallest eigenvalue of a sparse matrix 
of dimension i A+ ^Y~ 1 \ where A is the number of atoms in W lattice sites. 

2.2 Quantum computation with a ID optical lattice 

In order to use the Mott-insulator phase for building a quantum register, it 
is advantageous to use atoms with two degenerate ground states \g a ) and 
| <7f,) which are coupled to each other with a Raman transition via an excited 
state [10]. In order to trap both 'species' of atoms simultaneously, an optical 
lattice formed of two counterpropagating laser beams with parallel linear 
polarization vectors is needed. The result of this configuration is that the 
atoms are trapped in two overlapping optical lattices with polarizations cr+ 
and er_ each of which can trap one of the two atomic ground states. 

If we denote the creation operators of atoms in the states \g a ) and |<?b) by 
cr and b' , respectively, we can write the interaction Hamiltonian (7) as 



— <« - 1) + U ah n>\ + —n\(n\ - 1) 
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- ( J ? a la l+ i + J-b\b l+1 + Jf atSi + h.c.) . (11) 

i 

Here we encounter several couplings amongst atoms of the same species (col- 
lisional couplings U aa , Ubb and tunneling rates J a , J b ) as well as coupling 
of atoms of different species (collisional coupling U a b and effective Raman 
coupling J R ). 

We assume that initially all atoms are in the ground state \g a )- This is the 
state which from now on will be denoted by |0) whereas the second ground 
state \gt) will serve as the logical |1). 

Single- qubit rotations 

From the above definition of the logical states it is clear that single-qubit 
rotations can be performed by a Raman process. For example, starting at 
|0) and performing half a Raman cycle leaves us with an atom in state |1), 
whereas shorter interactions produce superposition states between |0) and 
|1). If one starts from an arbitrary superposition, the application of a full 
Raman cycle generates an effective <r x -rotation. In addition, a simple Rabi 
cycle of 2tt applied to the ground state \gb) is equivalent to changing the 
phase of the logical |1) by it and thus represents a a z operation. 

Two-qubit gates — controlled-phase gate and square root of the swap 
operator 

In order to realize two-qubit operations, we need to couple two atoms in neigh- 
boring lattice sites to each other in a controlled way. Let us assume that the 
lattice is sufficiently deep such that without external changes tunneling be- 
tween lattice sites is prohibited, hence J a and J b can be initially neglected. A 
tunneling coupling between neighboring sites can be switched on by applying 
an additional standing wave perpendicular to the lattice with a waist that 
should not exceed the size of two lattice sites (see Fig. 6). Depending on the 




Fig. 6. The tunneling interaction between neighboring lattice sites can be switched 
on by lowering the potential barrier between the sites. 
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circular polarization (er + or er_) of the additional laser beam the tunneling 
strength for both atomic species can be tuned independently. 

For the two-qubit operations to work, the tunneling coupling must al- 
ways be much smaller than the collisional coupling, even when the potential 
barrier has been lowered. Hence, we require U S> J for all times during the 
interaction process. This amounts to the fact that the system stays in the 
eigenspace of the collisional terms U aa and Ubb, hence in a number state. 
This eigenspace is in fact degenerate with respect to the occupation numbers 
n a ,b = and n a j — 1 since the collisional energy is obviously zero in both 
cases. There is an energy gap from this degenerate subspace to the state with 
two atoms of the same species per lattice site which means that those states 
can be adiabatically eliminated from the evolution. If, in addition, there is a 
large inter-species collisional coupling U a b, then the same energy gap persists 
even for states in which two atoms of different species occupy one lattice 
site. Thus, every state with more than one atom per lattice site, regardless 
of their species, can be adiabatically eliminated thereby leaving us with a 
degenerate eigenspace spanned by the logical states |0) = \n a = l,n& = 0) 
and |1) = \n a = 0, n& = 1). This constitutes our well-defined computational 
space. 

In what follows, we will use the notation introduced in [10] by denoting the 
state of the atomic population in two lattice sites by \n\, n\\r? a ^ n\). Then, for 
example, a state of two atoms in their respective ground states \gb) is given 
in this notation by |01; 01) and represents the logical two-qubit state |11). 
Suppose now we were to lower the potential barrier between two neighboring 
sites only for the atom in ground state \gb) which can be done by choosing an 
appropriate polarization of the incident laser beam. Then, according to the 
Hamiltonian (11), this state can couple to only two other states, 1 02; 00) and 
1 00; 02), with two atoms simultaneously at one lattice site. The Hamiltonian 
(11) can be written in the basis {|01; 01), |02; 00), |00; 02)} as 



This Hamiltonian effectively corresponds to a ^-system with ground state 
|01 ; 01) and excited states |02; 00), 1 00; 02) coupled by an effective Rabi fre- 
quency — J b /2 and detuned by Ubb (see left figure in Fig. 7). Assuming that 
the detuning is large, the system remains in the ground state and acquires 
a phase —2 J dt (J 6 ) 2 /Ubb- However, in the same manner the logical states 
|01) = 1 1, 0; 0, 1) and 1 10) = |0, 1; 1, 0) acquire phases. This is due to their in- 
teraction with the states |1, 1; 0, 0) and |0, 0; 1, 1), respectively. For example, 
in the basis {|1, 0; 0, 1), |1, 1; 0, 0)} the Hamiltonian (11) simply reads 




-J b -J b 
J b U hb 
J b Ubb 




(12) 




(13) 
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I2,0;0,0> I0,2;0,0> ll,l;0,0> I0,0;1,1> 




-J b 



I0,1;0,1> I0,1;1,0> ll,0;0,l> 



Fig. 7. Effective 3-level system coupling the state |0, 1;0, 1} to its excited states 
(left figure), and effective 4-level system coupling |0, 1; 1, 0) and |1, 0; 0, 1) to their 
respective excited states (right figure). 



which leads to a phase shift of J dt (J b ) 2 jU a b- The latter can be compensated 
for by applying a single-qubit rotation on the state |1) on both lattice sites. 
The overall effect is to introduce a phase shift 



T 

= 2 M 



(j b ) 2 (j b f 

U ab ~ 



(14) 



to the logical state |11) while keeping all other logical basis states unchanged. 
This is the controllcd-phase gate which, together with the single-qubit rota- 
tions, constitutes a universal set of operations for quantum computing. 

Similarly, we can choose to act upon the logical states |01) and 1 10) which 
contain one atom of each species in neighboring lattice sites. In the adiabatic 
approximation, the Hamiltonian connecting these two states is (in the basis 
{|1, 1; 0, 0), |1, 0; 0, 1), |0, 1; 1, 0), |0, 0; 1, 1)}, see right figure in Fig. 7) 



H a b 



( Uab -J a \ 

-J a -J b 
-J b -J a 

-J b -J a Uab 



(15) 



J 



with the solution that the following transformation on the logical states |01) 
and 1 10) is achieved: 



|01) i ^ e~ iv [cos/|01) -isin/|10)] 



1 10) 



-up 



-isin/|01) + cos 1 10)] 



with 



9 : 



I 

I 



b\2 



(It 



(j a r + (j ) 

Uab 



1 

'J 



dt 



J a J b 



(16) 
(17) 



(18) 
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The effective Hamiltonian is therefore just 

ff eff = -/|(|10}(01| + |01)(10|), (19) 

which, for / = 7r/4 is also known also the square root of the swap operator. 
Additionally, the logical states 1 00) and |11) acquire phases —2 J dt (J a ) 2 /U aa 
and —2 J dt (J b ) 2 /Ubb, respectively, as discussed before. Thus, the two qubits 
pick up only an overall phase if one chooses (J a ) 2 /U aa + (J b ) 2 /Ubb = [{J a ) 2 + 
(J b ) 2 ]/U a b- This overall phase can later be reversed by applying appropriate 
single-qubit rotations on both lattices sites. 

These examples show that it should be straightforward in principle to 
implement different types of single- and two-qubit operations with very few 
laser pulses (typically just one). The drawback is that for all above consider- 
ations the adiabaticity condition must be fulfiled. Thus, the couplings J a and 

are assumed to be small compared to the collisional couplings U aa , Ubb, 
and U a b which results in unwanted long gate evolution times. For example, in 
Ref. [10] it has been estimated that, for currently measured collisional cou- 
plings of O(lkHz) [11], the gate operation time for nontrivial gates such as 
the above-described controllcd-phase gate with an error rate less than 10~ 3 
is roughly 100ms. This time scale is far too long to render quantum com- 
putation useful with this scheme. A way to circumvent this problem and to 
drastically reduce gate operation times is to relax the adiabaticity condition 
and to note that even without adiabatic evolution there are certain instances 
in which the atomic population returns completely to the logical space in 
which we have started. The drawback here is that laser amplitudes and pulse 
durations have to be stabilized much more precisely than in the adiabatic 
regime. Despite that, it seems that there is much potential in this and other 
proposals for quantum computing on an optical lattice. 

Moreover, this time scale is already of the order of the currently possible 
trapping lifetime of atoms in recent experiments (see Sec. 4 for a detailed 
discussion about trapping losses). Hence, only a few gate operations can be 
performed before the atoms are lost from the optical lattice. 

2.3 Experimental realization with atom chips 

A particularly interesting way of implementing the above ideas of atomic 
registers is by using atom chips (for reviews on this exciting subject, see for 
example [12, 13]). The basic idea here is to trap cold atoms in one of their 
low-field seeking hyperfine ground states in the combined magnetic fields 
of a current-carrying wire and a constant transverse bias (see Fig. 8). The 
radial magnetic field of the wire is superimposed on a constant homogeneous 
field from the side thereby creating a line of zero magnetic field parallel to the 
wire. The atoms are thus trapped in a tubular region whose distance from the 
wire can be tuned by adjusting the relative strengths of the two overlapping 
magnetic fields. However, as depicted in the inset of Fig. 8, the atoms can flip 
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Fig. 8. Principle of the magnetic guide for atoms (see text for details). 

their magnetic sublevels once they reach the region of zero magnetic field, 
hence yet another bias field has to be applied, this time parallel to the wire, 
to prevent the atoms from doing that. 

Once the atoms are trapped and evaporatively cooled to form a Bosc- 
Einstcin condensate, the atom cloud has typically a length of 100/mi and a 
diameter of 2/mi, hence it forms a long and thin cigar-shaped object which 
can be treated as quasi-lD. The idea is now to confine the cloud between 
two highly reflecting mirrors that form a microcavity in the standing wave of 
which the atoms experience the periodic light potential we were envisaging 
earlier (Fig. 9). The idea is then to move a string of atoms [14] in a condensate 




Fig. 9. Proposed scheme for combining atom chips with optical cavities for gen- 
erating one-dimensional optical lattices. 
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state into the cavity region to perform the transition to the Mott-insulator. 
The current experimental status is that experiments have been performed 
in which a cloud of trapped atoms (not a BEC yet) has been moved in a 
'conveyor- belt' fashion across the surface of a miniaturized atom chip. The 
next step will be to integrate two adjacent optical fibres with polished ends 
onto the atom chip to create the standing-wave potential. 

An obvious performance limitation is given by the time the atoms actu- 
ally spend in their respective trapped states above the wire. Several noise 
sources, both technical and fundamental, can cause spin flips from a trapped 
to an anti-trapped magnetic sublevel that will prevent the atoms from staying 
indefinitely above the wire surface. The influence of the predominant noise 
source — magnetic fluctuations caused by absorption in the current-carrying 
wire — will be investigated in more detail in Sec. 4. 

3 Photonic realisation — passive linear optics and 
projective measurements 

In a rather different setting compared to Sec. 2 we can consider photons as 
the carriers of quantum information. Photons constitute an alternative to 
atoms — being qubits at rest — as they are massless particles and therefore 
move at the speed of light which earned them the name 'flying' qubits. They 
will eventually be part of larger networks and are considered to be vital 
in transporting quantum information over longer distances. Moreover, it is 
believed that photons could even be used to perform quantum operations 
themselves. It is therefore vital to know what kind of operations can be done 
with photons and how they are implemented. 

3.1 Qubit encoding and single-qubit operations 

First of all, it is necessary to define which photonic degrees of freedom wc 
would like to choose to encode our qubits in. One possibility would be to en- 
code the information in the polarization state (horizontal or vertical) [15, 16], 
another one the superposition states of one photonic excitation in two modes 
[17] . Another, seemingly complementary encoding would use the photon num- 
ber or Fock states of a photon. All the mentioned possibilities have their 
advantages and disadvantages. For example, single-qubit rotations are triv- 
ially implemented in the polarization basis because they are just performed 
by A/4- or A/2-plates. However, two-qubit operations such as the controlled- 
phase gate are impossible to implement in this way using wave plates. 

In what follows we will look at an alternative encoding in which the qubits 
are defined by photon numbers. That is, the logical will be the vacuum 
state of the electromagnetic field, and the logical 1 will be a single-photon 
Fock state. Then, let us consider a simple example of a two-qubit gate, the 
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controlled-phase gate defined by the truth table 2. That is, only the basis 
state containing one photon in each mode will pick up a phase tp, all other 
basis states are left unchanged. 

Here we immediately encounter a problem in that, in order to realize this 
quantum gate, two single photons have to interact with each other sufficiently 
strongly to produce the desired phase shift. In fact, nature is not so kind 
as to allow us to do this easily. Consider for a moment standard quantum 
electrodynamics. The lowest-order Feynman diagram that contains a photon- 
photon interaction is depicted in Fig. 10. The interaction strength is of fourth 
order in the fine structure constant a w 1/137 and therefore negligible. 




Fig. 10. Lowest-order Feynman diagram leading to a photon-photon interaction. 

A seemingly promising alternative is provided by nonlinear materials ex- 
hibiting a or Kerr nonlinearity. The constituents of those materials re- 
spond nonlinearly to an external electromagnetic field, and after tracing out 
the matter degrees of freedom leaves behind an effective nonlinear interaction 
between photons. However, these natural nonlinearities are still too small to 
be of any use since only phase shifts of the order of 10~ 8 can be achieved. 

3.2 Measurement-induced nonlinearities 

A possible way out of this dilemma can be found in the use of so-called 
measurement-induced nonlinearities. The idea behind it is perfectly simple. 
Suppose we wanted to act on a pure single-mode state \ip) with a nonlin- 
ear operator. In order to see this, we mix our signal state with another pure 
single-photon state at a beam splitter and perform a suitable projection mea- 
surement at one output port of the beam splitter. The result will be a con- 
ditional nonlinear operator acting on the signal mode. For example, let us 
consider the situation depicted in Fig. 11 in which a single-photon Fock state 
acts as our auxiliary (or ancilla) state. It it well-known that a beam splitter 
acts as an SU(2) group element on the photonic amplitude operators [18]. 
Hence, the operator associated with the action of the beam splitter can be 
written in the two-mode representation of the Pauli operators (see Sec. 1) as 
[19] 
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Fig. 11. Single beam splitter acting conditionally as a nonlinear operator. 



U12 = T" 1 e~- R * a 2 ai e i ^ a2 !T~ fl2 (20) 

where 1,2 label the modes impinging on the beam splitter. Consider now the 
situation in which the auxiliary mode is prepared in a single-mode Fock state 
and a single photon is detected in the same output port. Then the resulting 
conditional operator acting on the signal mode reads 

% = (1 2 |L> 12 |1 2 ) = T^- 1 [\T\ 2 - jhIRI 2 } . (21) 

It is clear from the expression (21) that this operator represents a nonlinear 
evolution since it is inherently quartic in the photonic amplitude operators 
a\ . This is an effective nonlinearity since the overall evolution, counting 
the contributions from all possible measurement outcomes, is still perfectly 
linear. One observes furthermore that the operator Y\ is independent of the 
signal state, in fact, the nonlinearity is created by the measurement process 
only. Therefore, following the notation commonly used in quantum optics, 
this represents quantum-gate engineering — as opposed to quantum-state 
engineering — or, in other words, quantum-state engineering of arbitrary 
states. Note also that, on the other hand, in general the success probability 
p = H^ilV') II does depend on the chosen signal state \ip). However, if the beam 
splitter parameters are chosen such that Y^Yi\ip) cx |V>), i-e. the conditional 
operator is proportional to a unitary operator (hence a quantum gate), then 
the success probability is state-independent as well. This important fact will 
be used later on to actually construct quantum gates with maximal success 
probability. 

In order to go one step further in the development of conditional quantum 
gates for photons, we have to look at larger beam splitter networks consum- 
ing more than just one auxiliary or ancilla state. To do so, we note that every 
V(N) transformation of N photonic amplitude operators can always be real- 
ized by a triangular-shaped network of at most N(N — l)/2 beam splitters 
and some phase shifters [20]. Hence, we need to generalize our arguments 
to whole networks of beam splitters. There seem to be essentially two ways 
of analyzing U(JV) networks. One is to look at the unitary operator U\2...n 
acting on N modes in terms of its Euler decomposition analogous to Eq. (20). 
Although this is not completely impossible, it is very tedious indeed. More- 
over, the Euler decomposition has only been calculated for N = 3, 4 [21]. An 
alternative way has been developed in [22] and uses well-known techniques 
from bosonic operator algebras. 
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For this purpose, we write the input state — still assumed to be a single- 
mode state, the theory is analogous for multi-mode states — in a functional 
form as 

|V)=/> t i)|0)=^^=(alr|0), (22) 



where the c m are constrained in such a way that |c m | 2 = 1. Analogously, 
we write the auxiliary state |A) and the state |P), respectively, in product 
form as 

|A)=T]W lo)^" 1 , |P) = rf-4=|0)® JV - 1 . (23) 



t =2 f = \ V^7- 

Here nij is a non-negative integer that represents the number of photons 
initially in mode i, and rij is the number of photons in the projected mode 
j. The beam splitter network is represented by a unitary A x A-matrix A, 
under the action of which the amplitude operators transform as 

a^A+a, k ] ^A T kK (24) 

Combining all these definitions we derive the (un-normalised) output state 
after mixing at the beam splitter network and projecting onto |P) as 

IV') oc (P|f/i2...iv|A) ® |V) 

N /- x„. / N \ mi /N \ 

= N - i9 <o\ n -W=f (E^ a t / e^ k"-^) 

We can see immediately from Eq. (25) that the effect of the beam splitter 
network is to mix the photonic creation operators of signal and auxiliary 
modes. At this point we use a well-known ordering formula 

a,F(a,tf)]=-^F(a,tf) (26) 
to rewrite Eq. (25) in the convenient form 

W) cc *-*<0| n / (X>«l ) \0r N . (27) 

i,j=2 V m *- n 3- \fe=l / / 

In the following we quote some basic results that follow immediately from 
Eq. (27). Suppose all A — 1 auxiliary modes are prepared in single- photon 
Fock states, and all A — 1 photodetectors find only the vacuum state at the 
respective output ports. Then the output is proportional to the A — 1-fold 
application of the creation operator (aj) w_1 . Similarly, if all auxiliary modes 
are prepared in the vacuum state and all photodetectors find exactly one 
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photon each, then the output is proportional to a^ -1 . These results should 
not be surprising if one remembers that we are merely selecting a particular 
measurement result from a passive linear operation which preserves photon 
numbers. Somewhat more interesting is the situation in which all auxiliary 
modes are occupied by single photons and all detectors register a photon each. 
Then the associated conditional operator is a polynomial of degree N — 1 in 
the number operator hi, P/v_i(ni). The last result is particularly important 
since it allows us to act upon the (unknown) coefficients of an iV-dimensional 
signal mode independently. 

3.3 Construction of simple quantum gates 

"The permanent doesn't really interact well with linear algebra." 
— Mark Jerrum, private communication 

The theory outlined above presents the framework which we will use to 
construct particularly simple and useful quantum gates. For simplicity, we 
first restrict ourselves to single-mode gates and give an outlook to multi- 
mode gates later. We also assume that we want to operate within Fock layers, 
i.e. we merely change phases of some expansion coefficients. The simplest 
nontrivial operations can be generated with signal modes that contain up to 
two photons. 6 Let us therefore consider the nonlinear phase shift gate C v 
which is defined by its action on a three-dimensional single-mode state as 

c |0) + ci|l) + c 2 |2) % co|0) + ci]l> + e^c 2 |2) . (28) 

According to what we have said earlier, in order to be able to act upon 
the coefficients of this three-dimensional signal state independently, we need 
a second-order polynomial in the number operator. Hence, we immediately 
find a network that performs the sought gate operation. The result is shown 
in Fig. 12. Note the triangular-shaped arrangement of the beam splitters as 



* lfo ul > 



Fig. 12. SU(3) network realizing the nonlinear phase shift, 
predicted in [20]. 



6 Note that all operations of the type co|0) + Ci|l) t— > c o |0) + e' v Ci\l) can be 
realized deterministically, since they are merely phase shifts. 
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Let us now analyze this network along the lines described above. We 
denote by A the unitary (3x3)-matrix associated with the SU(3) network. 
The conditional operator Y acting on the signal mode is then 

= coper A(1|1)|0) + ciper A|l) 
+c 2 (2/l llP er A - yl? lP er A(l|l) + 2A 12 A 21 A 13 A 31 )\2) (29) 

where per A denotes the permanent of the matrix A and per A(l|l) its prin- 
cipal subpermanent. In order to realize a nonlinear phase shift, we need to 
fulfil the relations 

perA(l|l) = perA, per A(l|l) [e iv + A\ x - 2A U ] = 2A 12 A 21 A 13 A 31 . 

(30) 

It turns out that there are infinitely many solutions to these equations which 
differ in their respective success probabilities. We have found numerically 
that the maximal probability is 1/4 [22] . 

Permanents 

The appearance of permanents in these problems is generic. The permanent 
of an (n x n) matrix A is a generalized matrix function, defined as 

n 

per A- tl^i ( 31 ) 

{a,}£S„ i=l 

where S n is the group of cyclic permutations [23]. Permanents naturally ap- 
pear in combinatorical problems, graph theory and related subjects. In our 
context they 'count' the ways of redistributing N single photons through an 
SU(7V) network to yield exactly N single photons at the outputs, i.e. 

perA = n ®(1\Ui 2 ...n\1)® N ■ (32) 

As one can see from Eq. (29), the unitary network has to be adjusted such 
that the permanent and certain subpermancnts of the associated unitary 
matrix fulfil certain relations. Furthermore, the overall success probability is 
just |per A(l|l)| 2 . Therefore, the optimal network is obtained by maximizing 
some permanent under a number of given constraints. 

This is, in fact, a nontrivial task which is partly due to the fact that 
the algebraic property of a matrix being unitary does not imply any ma- 
jor simplifications for computing permanents. 7 One of the few known fact 
is that permanents of unitary matrices are bounded from above. This is a 
consequence of the Marcus-Newman theorem [23] which states that for all 
(m x n)-matrices U and (n x m)-matrices V the inequality 



7 This is in stark contrast to the determinant which is just the product of the 
eigenvalues of a matrix. 
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|per UV| 2 < per UlTper VV* (33) 



holds which reduces, when setting V = I and regarding U as being unitary, 
to |perU| < 1. This is, of course, what one suspects if the permanent is 
supposed to be related to a success probability. 

Another interesting fact to note is that computing permanents of ma- 
trices of increasingly larger size is a computationally hard problem in the 
sense that the computing time scales exponentially with the size of the ma- 
trix. Computer scientists say this problem is NP-complcte. It seems that this 
computational problem is related to our quest to design quantum networks 
that would eventually be able to solve problems in polynomial rather than 
exponential time. That is, we have to solve an exponentially hard problem 
first by designing networks in order to avoid it later when applying them 
to unknown quantum states in the course of a quantum computation. Note 
however, that this argument strictly applies only if one is interested in the op- 
timal network. For most practical purposes it is sufficient to approximate the 
permanent which can be done efficiently with a quadratic number of steps. 

Dimension of the auxiliary state 

What we have said so far about the nonlinear phase shift gate regarded the ex- 
plicit construction of conditional nonlinear operators. In turns out, however, 
that this way of designing quantum gates, although being algorithmically 
transparent, does not yield optimal networks in terms of resources. For ex- 
ample, there is an alternative network proposed in Ref. [24] that also realizes 
the nonlinear phase shift but with one fewer beam splitter (see Fig. 13). As 
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Fig. 13. Nonlinear sign shift gate with two beam splitters. 



we will see later, the best networks for multi-dimensional quantum gates in 
terms of their respective success probabilities are not derived in this manner. 

What seems to influence the required resources most is the dimensionality 
of the used auxiliary state. Comparing Fig. 12 and Fig. 13 one realizes that 
the Hilbert spaces spanned by the auxiliary modes are different in both cases. 
Starting from a product state of the form |11) as in Fig. 12, by beam split- 
tings we actually span a three-dimensional space formed of the basis states 
{|20), 1 11), 1 02)}. On the contrary, with an initial product state of the form 
1 10) as in Fig. 13 we merely span a two-dimensional space with basis states 
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{|10), |01)}. From this observation we can conclude that it is the dimension- 
ality of the Hilbert space of the auxiliary modes that determines our ability 
to design quantum networks with the least possible resources. 

3.4 Multi-mode gates 

Similarly to the construction of single-mode quantum gates such as the non- 
linear phase shift, we can proceed to more complex networks by referring 
to the results obtained before. The idea is rather simple. Let us restrict our 
atttention again to quantum operations within a certain Fock layer, i.e. a 
subspace of the total Hilbert space of a multimode system with fixed photon 
number. Then we construct an M-mode quantum gate in the following way: 

1. first feed the M modes into a generalized Mach-Zehndcr interferometer 
with M input and output ports (2M-port for short), 

2. then act upon each output mode with an appropriate single- mode gate, 

3. and finally rccombine the M modes at another 2M-port. 

Note that the first and last steps are done deterministically. 

Consider for example the controlled-phase gate whose truth table is given 
in table 2. Because it represents a two-mode gate acting in the Fock layer 
with total photon number equal to two, this gate serves as the prime example 
for our statement. In terms of photonic amplitude operators, the conditional 
operator associated with the controlled-phase gate is given by 

C v = 1 - (1 - e^)n x ni , (34) 

where we have already written the operator in the simplest form that matches 
the dimensionality of the signal-mode Hilbert spaces. 8 If we decompose the 
operator C v into a product of two unitary operators associated with the 
Mach Zchndcr interferometer sandwiching some other nonlinear operator we 
find that these nonlinear operator indeed forms a tensor product N± ® N2 of 
two single-mode conditional operators corresponding to nonlinear phase shift 
each of which has the form 

= 1-^(1-^)^(^-1), i = 1,2. (35) 

Thus, combining all results we have obtained so far, we end up with a net- 
work depicted in Fig. 14. With the knowledge about the maximal success 
probability for a nonlinear sign flip (a nonlinear phase shift with ip — tt) 
being 1/4, the corresponding success probability of the controlled-phase gate 
will be 1/16. 

As we have already mentioned, this is not the highest probability one 
can actually achieve with linear optics and photodetection. The 'record' up 

8 Here we assume that the two-mode signal state is indeed only formed by the 
basis states { |00) , |10), |01), |11}} without contributions from multi-photon states. 
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Fig. 14. Controlled-phase gate as two nonlinear phase shifts inside a balanced 
Mach-Zehnder interferometer. 



until now is set by Ref. [25] where a network has been found that works 
with a probability of 2/27. This network does not have the structure we have 
discussed here. 



3.5 Conditional dynamics and scaling of success probabilities 

In connection with the original proposal on using post-selection of sub- 
dynamics it has been suggested that the use of conditional dynamics could 
eventually yield success probabilities that are arbitrary close to unity. What 
is meant by conditional dynamics is that depending on a certain measure- 
ment outcome (or pattern for more detectors) the resulting output state is 
either retained (in case the measurement gave the wanted answer) or post- 
processed using the knowledge about the transformation due to the detection 
result. We will give a simple argument that within the framework of qubits 
being encoded in photon numbers this cannot hold. 

For this purpose, let us return to the single-mode gate described earlier. 
Suppose therefore we had a single-mode quantum state of the form co|0) + 
Ci|l)+C2|2) at hand and we send it through the arrangement of beam splitters 
as depicted in Fig. 13. After the first beam splitter which is fed with an 
auxiliary single photon the output can contain as much as three photons 
in one output port. Hence, there is a certain probability that two or three 
photons are found in the photodetector. Let us look at the extreme case in 
which all three photons end up in the detector. By conservation of photon 
numbers, the output state contains only vacuum. But that means that we 
have lost all information about the signal state which was originally encoded 
in the weights a of the superposition. If we detected two photons instead 
we would lose the information contained in cq. Thus, we can conclude that 
conditional dynamics breaks down if the total of number of photons measured 
in all photodetectors exceeds the total number of photons in the auxiliary 
state. Having said that, the obvious way of rescuing this situation is by noting 
that we have taken the measurements to be projective measurements. If in 
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future one finds a way of performing non-destructive (QND) measurements, 
one has immediately circumvented the information loss since all that has 
happened is a unitary transformation of the signal state. 

Up until now, we have to live with projective measurements and it makes 
sense to ask what probabilities can be achieved in principle for certain gates 
and how they scale with increasing dimensionality of the signal state. We have 
already seen that the nonlinear sign flip can be realized with p = 1/4. In view 
of the algorithm to design multi-mode gates using generalized Mach-Zehndcr 
interferometers, we could ask what the respective success probabilities are 
if we started from an N + 1-dimensional signal state of the form Co|0) + 
Ci 1 1) + ... + cn\N) and we would like to design a network that realizes a 
sign flip on the coefficient cat. Then in turns out that if one used an TV- 
dimensional auxiliary state with the lowest possible photon numbers, i.e. 
photon numbers in the range . . . N — 1, there are networks that achieve a 
the desired transformation with a probability of exactly 1 /N 2 [26] . 

4 Decoherence mechanisms — QED in causal dielectric 
media 

"Please mind the gap between the hope and the reality." 
- Smoke #2 

4.1 Decoherence mechanisms affecting atoms and photons 

Trapping losses near current- carrying wires 

As discussed in Sec. 2, the time the trapped atoms spend above the current- 
carrying wire is limited by several loss and heating mechanisms. In Ref. [13] 
one can find a table with several loss mechanisms relevant for atom chip 
experiments. We summarize the most important ones in table 3 and give 
their approximate associated lifetimes. Most of the loss mechanisms are in 

Table 3. Loss mechanisms for trapped atoms above current-carrying wires 
mechanism lifetime 

thermally induced spin flips 1 — 10s 
technical noise > 10s 

background collisions > 10s 



fact avoidable as they are merely technical noise or account for experimental 
imperfections. Technical noise can cause both heating and spin flips. Heating 
processes are associated with low-frequency noise that causes the atoms in 
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the trap to be excited into higher vibrational modes. Spin flips, on the other 
hand, are induced by rf noise. Suppose an atom stayed Is over the wire. The 
Rabi frequency that induces the spin flips is therefore 0.5Hz, which in turn 
is approximately given by 10 6 Hz/G. Therefore, the rf magnetic field must be 
less than 1/iG which, for an atom being 1mm away from the wire, amounts 
to current fluctuations of less than lfiA. For a current of 1 — 10A needed 
to produce the dc magnetic field this means that the power supply must be 
unusually stable. This challenging but achievable. 

However, some loss mechanisms are inherent in the way the atoms are 
trapped and can hardly be avoided. Most notable of these unavoidable losses 
is that caused by thermally induced spin flips which incidentally already gives 
rise to the biggest loss rate in current experiments [27]. In Figure 15 we show 
the five magnetic sublevels of the |.F = 2)-state of 87 Rb. In the experiment, the 
atoms are optically pumped into the |_F = 2,to = 2) sublevel in which they are 
trapped. Spin flips induced by fluctutating magnetic fields cause the atom to 
jump into lower-lying magnetic sublevels in which they are either lost due to 
gravitational forces (|F = 2,m = 0)) or expelled from the trap (\F = 2,m = — 1) 
and \F = 2, m = — 2}). These spin- flip transitions are caused by magnetic- field 
fluctuations that are unavoidable as soon as metallic or dielectric materials 
are located close to the atoms. This is quite obviously the case in atom chips. 
In order to see where these fluctuations come from, we have to look deeper 
into the statistical implications of macroscopic (quantum) electrodynamics. 
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IF=2,m=-l> (anti-trapped) 
IF=2,m=-2> (anti-trapped) 



Fig. 15. Zeeman splitting of the \F = 2) state of 87 Rb and some of the possible 
spin-flip transitions. 



Gate fidelity with imperfect linear optical elements 

In Sec. 3 on all-optical realisations of quantum information processing with 
passive linear optical elements we have seen that the fidelity of a gate op- 
eration will depend crucially on the ability to generate single photons with 
very high efficiency and the amount of losses that are introduced by imper- 
fect passive elements and photodetectors. It is not difficult to imagine that 
the purity of the ancilla photons used to operate a quantum circuit will have 
enormous influence on the way the quantum gates operates. For example, 
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an ancilla state consisting of a mixture of vacuum and a single-photon Fock 
state actually corresponds to realising two different quantum gates at once or, 
in other words, a statistical mixture of their respective operations. Another 
example of a decoherence process is provided by the inherent absorption in 
the passive linear optical elements used to build up the quantum circuit. An 
absorbing beam splitter does not act unitarily on the electromagnetic field 
modes impinging on it, but rather as a general CP map, i.e. as a statistical 
mixture of the desired and some unwanted operations. 

All the loss and decoherence mechanisms sketched above can be described 
by a single and powerful theory which we will outline in this section. 

4.2 Field quantisation in causal media 

"Die ganzen Jahre bewufiter Griibelei haben mich der Antwort der 
Frage 'Was sind Lichtquanten' nicht naher gebracht. Heute glaubt 
zwar jeder Lump, er wisse es, aber er tauscht sich." 
— Albert Einstein, from a letter to Michele Besso (1951) 

Quantum electrodynamics (or QED for short) is a well-established theory, 
and one might wonder how it is possible to obtain fundamentally new results 
apart from those that are already known. The answer comes about when 
considering phenomenological electrodynamics and when trying to quantize 
it. 

Classical electrodynamics with media 

Maxwell's equations of the electromagnetic field in free space, i.e. without 
external sources, in temporal Fourier space can be cast in the following form: 



These equations have to be supplemented by appropriate constitutive re- 
lations which connect the electric field E(r,w) and the magnetic induction 
B(r, w) with the displacement field D(r, w) and the magnetic field H(r, w) 
which also carry information about the material. Assuming a purely dielectric 
medium, i.e. disregarding magnetic properties, one usually defines a polari- 
sation field P(r, ui) via 



V-B(r,w) 
V-D(r,a/) 

V x E(r,w) 

V x H(r,w) 



0, 
0, 

iwB(r, u) , 
— iwD(r, u) . 



(36) 
(37) 
(38) 
(39) 



D(r, uo) = e E(r, u) + P(r, J) . 



(40) 



In real space and assuming locally and linearly responding media, the polar- 
isation can be written as a temporal convolution 
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oo 

P(r, t)=e J dr X (r, r)E(r, t - r) + P w (r, t) (41) 
o 

with the linear susceptibility x( r , T ), the Fourier transform of which is related 
to the dielectric permittivity e(r, w) by 

e{T,w) = X {r,u) + l. (42) 

The permittivity is a complex function of frequency e(r, u>) = en(r,uj) 
-Hei(r, ui). The real and imaginary parts, which are responsible for dispersion 
and absorption, respectively are related to each other by the Kramers-Kronig 
relations 

e R {r lU )-l = ^V (43) 
i j u) — U) 

e I (r,.) = - 1 -V fd ^^- 1 (44) 

7T J Lj' — LU 

with V denoting the principal value. The complex permittivity is an ana- 
lytic function in the upper complex half-plane without zeros and satifies the 
relation 

e(r,-w*) = e*(r,a;). (45) 

Furthermore, it approaches unity in the high-frequency limit. These analyt- 
icity properties, which carry over to the Green function to be defined later, 
will be important for proving equal-time commutation relations (ETCR for 
short) between electromagnetic field operators. This list of properties can be 
shown to follow from causality. 

Equation (41) contains an additional, non-convolutive term Pjv(r, t). In 
fact, the existence of this term is crucial for the functioning of the following 
quantisation procedure. It serves as a Langevin noise source which is needed 
to preserve the ETCR. From statistical physics it is known that dissipation 
processes (as described by the imaginary part of the dielectric permittivity) 
are always accompanied by additional fluctuations. In order to see this more 
clearly, assume a damped harmonic oscillator with the solution for the ex- 
pectation values (a(t)) = (a(i'))e _r (* _t \ This relation cannot be valid in 
operator form since with increasing time the ETCR between a(t) and a^(t) 
would decay exponentially and eventually violates Heisenberg's uncertainty 
relations. Certainly, it is possible to add a Langevin force with vanishing ex- 
pectation value to the oscillator which takes care of the ETCR. The noise po- 
larisation PAr(r, t) is exactly such a Langevin force. Its strength is determined 
by the linear fluctuation-dissipation theorem [28]. Assuming that Pjv(r, t) is 
proportional to a Gaussian random variable f(r, u>) which is always possible in 
linear reponse theory, and noting that the linear fluctuation-dissipation theo- 
rem states that the correlation function (P(r, u), P(r', a/) is proportional to 
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the imaginary part of the response function — hence the dielectric permittiv- 
ity — we immediately deduce that the correct form of the noise polarisation 
must be [29, 30, 31] 



Pjv(r,i)=iy^£/(r, W )f(r,u,). (46) 

As a matter of fact, the f(r, u) plays the role of the fundamental variable in 
terms of which all relevant quantities can be expressed. Consider for example 
the Helmholtz equation for the electric field which is obtained by substituting 
Faraday's law (38) into Ampere's law (39) and using the constitutive relation 
(40) as 

LO 2 

V x V x E(r,w) -e{r, w)E(r, u) = /i w 2 PAr(r, ui) . (47) 

The solution to Eq. (47) is easily found to be 

E(r, lo) = ^lo 2 j d 3 s G(r, s, to) ■ Pat(s, to) (48) 
where the (tensor-valued) Green function satifies 

V x V x G(r, s, u) - ^2"£(r, w)G(r, s, lo) = S(r - s)U , (49) 

where U denotes the unit dyad. Equation (48) shows explictly how all field 
quantities can eventually be expressed, via the noise polarisation, in terms of 
the fundamental field f(r,w). As mentioned earlier, the dyadic Green func- 
tion G(r, s,cj), being a causal response function itself, has all the analytic 
properties the dielectric permittivity has, i.e. it is holomorphic in the upper 
complex frequency half-plane and satifies the relation 

G(t,s,-(J*) = G*(t,s,u). (50) 

Moreover, it fulfils a generalized Onsager-Machlup reciprocity [28] relation 
of the form 

G(s,r,w) = G T (r,s,w) . (51) 

From the partial differential equation (49) it also follows that the dyadic 
Green function satisfies an important integral relation, 



/ 



LO 2 

d 3 s -^ei(s, w)G(r, s, w) ■ G+(r', s, w) = ImG(r, r', w) , (52) 



which expresses the fluctuation-dissipation theorem. 
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Quantized electrodynamics with media 

The theory described above can now be quantized by regarding the funda- 
mental field f(r, u) as a bosonic vector field f(r, w) with the commutation 
rule 

f(r, w), f V, u/)l = S(r - y')5{uj - J)U . (53) 

Then, for example the Fourier component of the electric field operator can 
be written, on using Eq. (48), as 

E(r, oj) = J d3s A/e/(s,w)G(r, s, u) • f(s, uo) . (54) 

The corresponding Schrodinger operator is obtained by integrating over all 
frequency, 

oo 

E(r) = J cL;E(r,w) + h.c. . (55) 
o 

The so-defined operators fulfil all the requirements one imposes on a sta- 
tistical quantum field theory. Firstly, the ETCR between electric field and 
magnetic induction can be shown to be [29, 30, 31] 

E(r), B(r')l - -— V x 6(r - r')U . (56) 
J e 

Secondly, the theory is — by construction — also consistent with the 
fluctuation-dissipation theorem which we can cast in the form 

(0|E(r,w)E t (r',w / )|0) = — ^jImG(r, r', u)5(u - w') . (57) 

7T£qC 

Equation (57) tells us that the (vacuum) fluctuations of the electric field 
are given by the imaginary part of the Green function as one would expect. 
Finally, the Hamiltonian describing the electromagnetic field in the presence 
of absorbing matter is given by 

oo 

H = j d 3 r J du)fuvf f {r,u))-i(r,oj) (58) 
o 

which is diagonal in the fundamental bosonic vector field f(r, w). Maxwell's 
equations then follow from the Heisenberg equations of motion of the above- 
defined electromagnetic field operators. 

Extensions of this theory to spatially anisotropic materials (for which the 
dielectric permittivity is a symmetric tensor) and media with gain can be 
found in [30, 31]. Magnetic materials, including the important class of mag- 
netodielectrics or left-handed media, can be treated in a completely analogous 
way [32]. We will, however, not elaborate more on this subject as it is not 
needed for our present purposes. 
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4.3 Thermally induced spin flips near metallic wires 

The above-described theory is now perfectly suitable to describe magnetic- 
field fluctuations radiated by metallic structures which cause spin flips be- 
tween trapped and anti-trapped (or non-trapped) magnetic sublevels of 
atomic hyperfine ground states. The experiment in [27] used a coated wire 
to trap 87 Rb atoms in their \F = 2,m = 2) magnetic sublevcl. The atoms 
are located at a distance r from the surface of the wire with radius a 2 (see 
Fig. 16). 



layer 3 
(vacuum) 



x 




T 

Fig. 16. A two-level atom is located at a distance r from the surface of a coated 
wire that runs along the 2-direction. 

The interaction of a neutral atom positioned at some point ta with the 
electromagnetic field described by the Hamiltonian (58) via the Zeeman in- 
teraction Hz — — /i • 1S(ya) of the atom's magnetic moment fx with the 
magnetic field causes the spin of the atoms to flip. Since the atoms are in 
their respective (electronic) ground state, it is safe to assume that the an- 
gular momentum is zero. Moreover, the nuclear magnetic moment is small 
compared to the Bohr magneton \xb- This leaves us with an expression for 
the atom's magnetic moment associated with the transition \i) — > |/) as 
ft = fi\i)(f\ + h.c. where the magnitude \x is given by the (electronic) spin 
matrix element /x = 2^s(i|S|,f). If we define the spin flip operators a = 
satisfying the angular-momentum algebra [<t,<t z ] = a, we can write the total 
Hamiltonian in rotating-wave approximation as 
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oc 

H = j d 3 r J dufiu)f\r,u))-f(r,u;) + ^7iWfi<j z 
o 

00 r~h 

a^fA-Vx dui d 3 r^j 1/ £/(r, w)G(ta, r, u>) ■ f(r, ui) + h.c. 

J J C Z V 7T£ 



.(59) 



Using the Hamiltonian (59) in the Hciscnberg equations of motion for the 
spin flip operators, we find after some calculation that in the Markov ap- 
proximation the spin flip rate is given by [33] 



(/|S|i> • lm [V x V x G(v A , r A , uj fi )} ■ (i\S\f) . (60) 



Taking into account that the flip rate is thermally enhanced, i.e. not only 
spontaneous but also stimulated (induced) flips occur, the rate (60) has to 
be modified by a multiplicative factor ( e huJ fi/ kT — 1)- 1 which is proportional 
to T for high temperature. Furthermore, considering the initial state \i) to 
be the \F — 2, m = 2) magnetic sublevel and the final state |/) to be the state 
\F = 2, m= 1), an additional factor of 1/4 due to the expansion of the \F, re- 
states in terms of the eigenstates of electronic and nuclear spin operators is 
present. 

Now we are in a position to compare with experimental measurements. 
The experiment in [27] used a coated wire with a copper core of 185/jm radius 
and an aluminium coating with 55/zm thickness. The transition frequency 
between two neighboring magnetic sublevels was / = Ufi/2ir — 560kHz. The 
result of the calculations with our theory and the experimental data are shown 
in Fig. 17 where we have plotted the inverse spin flip rate, hence the average 
lifetime r of an atom in the trap. It is clear that our theoretical result will 
only serve as an upper bound on the lifetime since, as we have mentioned 
earlier, there are indeed other noise sources that influence the duration of 
stay of the atoms in their trap. Note, however, that the agreement is still 
astonishingly good which means that our earlier claim that most other noise 
source are negligible in this experiment, was indeed valid. 



4.4 Imperfect passive optical elements 

Another, seemingly unrelated, application of the field quantization in absorb- 
ing materials is found if we consider passive optical elements such as beam 
splitters under realistic conditions, i.e. including losses in the beam splitter 
material. As noted earlier, a lossless beam splitter acts upon the amplitude 
operators of the incoming waves as an SU(2) group element [18]. That is, two 
amplitude operators d\{oj) and a 2 (w) are transformed into new amplitude 
operators 61 (u>) and 62 (to) as 
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T (s) 20 








20 40 60 80 100 

r (/ini) 

Fig. 17. Lifetime r of the trapped atom as a function of the distance r from the 
surface of the current-carrying wire. The solid line corresponds to the calculated 
lifetime according to the theory presented in the text. The crosses represent exper- 
imental data points. The dotted line is the result for a slab-geometry according to 
[34]. 



where T(w) is a unitary matrix, i.e. T(w)T + (cj) = J. For lossy beam split- 
ters, this is certainly not true since some of the electromagnetic radiation 
impinging on the beam splitter will actually be absorbed. This process can 
be accounted for by introducing an absorption matrix A(w) in order to sat- 
isfy the energy conservation T(a>)T + (cj)+ A(w)A + (u>) = I. Both matrices 
actually follow from a decomposition of the associated Green function into 
contributions associated with travelling waves moving to the left and right 
[35] and are therefore given by the complex refractive-index profile of the 
beam splitter material and geometric parameters. 

However, it turns out that, associating the material excitations with some 
bosonic variables gi(to) and gi (w) which are special linear combinations of 
the fundamental bosonic field variables f(r, ui) restricted to one spatial dimen- 
sion [35], the beam splitter acts unitarily on the combined set of photonic 
and material amplitude operators [36]. In fact, a lossy beam splitter can 
be represented as an SU(4) group element. That is, the 'four-vector' ct{uj) 
consisting of incoming photonic amplitude operators and the material exci- 
tations transforms unitarily into the 'four-vector' 0(uj) of outgoing photonic 
amplitude operators and some new matter operators, i.e. 

/§(w) = A(u)a(u) , A(u) A+ (w) = I . (62) 

The unitary 4x4- matrix A(ui) can be decomposed into block form as 
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A[UJ) ~ { -S(^)C- 1 (a ) )T(co) C( W )S- 1 (o ) )T(co) J [M) 

where C(a>) = ^T(lu)T + (lo) and S(lo) = yj A(w)A + (w) are commuting posi- 
tive Hermitian 2x2-matrices. The associated unitary operator U that realizes 
the transformation = A(u))a(u) =Wa(u))U is then obtained as 



f 7 T ] 

E/ = cxp < — i / duj [cr (w)] > 



(64) 



where the matrix #(w) is defined by e ~ l *( w ) = A(w). The operator U and 
the matrix A(uj) can be used as well to transform quantum states. For this 
purpose, let us assume that the density operator of the quantum state before 
the transformation is an operator functional of a(w), i.e. g — g[a(uj), 0.^(01)]. 
Then, transformation of the quantum state and tracing over the degrees of 
freedom associated with the beam splitter leaves us with 



g' = Tr \ua(u)U\ Ua\u)tP } . (65) 



In quantum information theory, decoherence processes are frequently de- 
scribed by the Kraus decomposition of the quantum-state transformation 
(65). That is, we seek an operator decomposition of the form 

# = Y^WiQW}, ^W^WW. (66) 

i i 

For this purpose, let us assume that the beam splitter is in its ground state, 
thereby neglecting thermal excitations which at room temperature and op- 
tical frequencies is a safe assumption. Let us also restrict ourselves to quasi- 
monochromatic radiation in which case we can drop all frequency dependen- 
cies. Then, calculating the partial trace in Eq. (65) in the coherent-state basis 
as 

e' = ^J d 2 a 3 d 2 a 4 W a3>a JW^ a4 (67) 
where the continuous- index Kraus operators are defined by 

W a3iCti = (a 3 ,a 4 \U\0 3 ,0 4 ), (68) 

we eventually find that the Kraus operators of an absorbing beam splitter 
are (exp -i& T = T) [22] 

W a3 , ai = cxp (-i[a f ] T <?> T a) cxp (-a+SC _1 Ta) exp (-a+a/2) . (69) 

One checks easily that these operators become unitary operators in the limit 
of vanishing absorption, i.e. when T is unitary and therefore <&t Hermitian, 
and S vanishes. We also see that the term exp(— a + SC _1 Ta) is responsible 
for absorption as it is a function of the annihilation operators only. 
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Other error sources 

Apart from absorption in the dielectric material, there are several other 
sources of errors in an experimental implementation. So far, we have as- 
sumed that the auxiliary states are perfectly pure states, especially we have 
implicitly assumed that we are able to produce single photons on demand 
with very high efficiency This, however, is an oversimplification. All realistic 
sources produce as most a mixed state of the form (1 — p)|0)(0| with 
efficiency p. The best single-photon sources to date achieve a single-photon 
efficiency of at most 80%. This means that the first beam splitter in Fig. 13 
does not always act as a first-order polynomial in n, but sometimes subtracts 
a photon from the signal state. 

Another crucial error source are the detectors themselves. In our mathe- 
matical description we have assumed that the photodetectors are perfect and 
can therefore be described by a projection operator |n)(n|. However, since 
they show losses themselves, we have to replace it by 

|n)<n| ~ n(n) = £ V ) k - n \k)(k\ (70) 

k=n ^ ' 

which constitutes a positive operator-valued measure (a generalization of a 
projective measurement). It describes the effect of absorption in the pho- 
todetector itself. The number r\ is commonly called the detector efficiency. 
Typical values for avalanche photo diodes are r\ w 0.3. 

We can measure the effect of these error sources by introducing the gate 
fidelity F as the overlap between the wanted result and the achieved 
density matrix g, i.e. F = (ip'\g\ip')- Since the so-defined fidelity still depends 
on the chosen signal state — because not all signal states are effected in the 
same manner by the errors — we define an average fidelity F by averaging 
over all possible signal states. 

Figure 18 shows the effect of non-unit single-photon efficiency p (left fig- 
ure) and detector efficiency rj (right figure) on the average gate fidelity. One 
can see that even a small inefficiency in either the detection process or the 
preparation process of the auxiliary state results in a (roughly linear) drop in 
the gate fidelity. In a different setting, in which the quantum information is 
encoded in the superposition states of one photonic excitation in two modes, 
a similar behaviour has been observed [37]. This result has important conse- 
quences for the scalability of such gates since quantum error correction only 
works if the individual errors in the network components does not exceed 
w 10~ 3 (see for example [38] and references cited therein). In other words, in 
order to be able to correct for errors, both discussed inefficiencies can be at 
most of the same order as the desired error rate which currently seems to be 
experimentally impossible. 
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Fig. 18. Average gate fidelity with imperfect auxiliary photons (left figure) and 
imperfect detectors (right figure). 
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